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Abstract 



We find that Polyakov lines, computed in abelian-projected SU(2) lattice gauge theory 
in the confined phase, have finite expectation values for lines corresponding to two units 
of the abelian electric charge. This means that the abelian-projected lattice has at most 
Z 2 , rather than U(l), global symmetry. We also find a severe breakdown of the monopole 
dominance approximation, as well as positivity, in this charge-2 case. These results imply 
that the abelian-projected lattice is not adequately represented by a monopole Coulomb 
gas; the data is, however, consistent with a center vortex structure. Further evidence is 
provided, in lattice Monte Carlo simulations, for collimation of confining color-magnetic 
flux into vortices. 
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1 Introduction 



The center vortex theory || and the monopole/abelian-projection theory || are two leading 
contenders for the title of quark confinement mechanism. Both proposals have by now 
accumulated a fair amount of numerical support. To decide between them, it is important 
to pinpoint areas where the two theories make different, testable predictions. In this article 
we would like to report on some preliminary efforts in that direction. 

Much of the numerical work on the center vortex theory has focused on correlations 
between the location of center vortices, identified by the center projection method, and 
the values of the usual gauge-invariant Wilson loops (cf. ||, Q). In the abelian projection 
approach, on the other hand, Wilson loops are generally computed on abelian projected lat- 
tices, and this fact might seem to inhibit any direct comparison of the monopole and vortex 
theories. However, it has also been suggested in ref. || that a center vortex would appear, 
upon abelian projection, in the form of a monopole-antimonopole chain, as indicated very 
schematically in Fig. |l|. The idea is to consider, at fixed time, the vortex color-magnetic 
field in the vortex direction. In the absence of gauge fixing, the vortex field points in arbi- 
trary directions in color space, as shown in Fig. [2|. Upon fixing to maximal abelian gauge, 
the vortex field tends to line up, in color space, mainly in the ±<r 3 direction. But there 
are still going to be regions along the vortex tube where the field rotates from the +03 
to the — cr 3 direction in color space (Fig. Upon abelian projection, these regions show 
up as monopoles or antimonopoles, as illustrated in Fig. |j. If this picture is right, then 
the ±27r monopole flux is not distributed symmetrically on the abelian-projected lattice, 
as one might expect in a Coulomb gas. Rather, it will be collimated in units of ±7r along 
the vortex line. We have argued elsewhere f| that some sort of collimation of monopole 
magnetic fields into units of ±7r is likely to occur even in the D = 3 Georgi-Glashow model, 
albeit on a scale which increases exponentially with the mass of the W-boson. On these 
large scales, the ground state of the Georgi-Glashow model cannot be adequately repre- 
sented by the monopole Coulomb gas analyzed by Polyakov in ref. 0. The question we 
address here is whether such flux collimation also occurs on the abelian-projected lattice 
of D = 4 pure Yang-Mills theory. 

The test of flux collimation on abelian-projected lattices is in principle quite simple. 
Consider a very large abelian Wilson loop 



corresponding to q units of the electric charge. The expectation values are obtained on 
abelian-projected lattices, extracted in maximal abelian gauge.0 If q is an even number, 
then magnetic flux of magnitude ±7r through the Wilson loop will not affect the loop. 

Abelian-projected links in Yang-Mills theory are diagonal matrices of the form = 
diag [exp (iA^ ) , exp (- i ) ] . 





or abelian Polyakov line 
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Figure 1: Hypothetical collimation of monopole/antimonopole flux into center vortex tubes 
on the abelian-projected lattice. 



Flux collimation therefore implies that W q (C) has an asymptotic perimeter-law falloff if q 
is even. Likewise, Polyakov lines P q for even q are not affected, at long range, by collimated 
vortices of ±7r magnetic flux. In the confined phase, the prediction is that P q = only for 
odd-integer q. In contrast, we expect in a monopole Coulomb gas of the sort analyzed by 
Polyakov that W q (C) has an area-law falloff, and P q = 0, for all q. 

Numerical results for abelian-projected Wilson loops and Polyakov lines must be in- 
terpreted with some caution in the usual maximal abelian gauge, due to the absence of a 
transfer matrix in this gauge. Since positivity is not guaranteed, these expectation values 
need not relate directly to the energies of physical states. For this reason, we prefer to in- 
terpet the abelian observables in terms of the type of global symmetry, or type of "magnetic 
disorder" , present in the abelian-projected lattice, without making any direct reference to 
the potential between abelian charges, or the energies of isolated charges. Following ref . || , 
let us introduce the U(l) holonomy probability distribution on abelian-projected lattices 




(1.3) 



for Wilson loops, and 




(1.4) 



for Polyakov lines, where 



1 



oo 




E e 



(1.5) 



is the 5-function on the U(l) manifold. These distributions give us the probability density 
that a given abelian Wilson loop around curve C, or an abelian Polyakov line of length T, 



Figure 2: Vortex field strength before gauge fixing. The arrows indicate direction in color 
space. 



Figure 3: Vortex field strength after maximal abelian gauge fixing. Vortex field strength 
is mainly in the ±03 direction 



Figure 4: Vortex field after abelian projection. 

respectively, will be found to have the value g G U(l) in any thermalized, abelian-projected 
lattice. The lattice has global U(l) symmetry (or, in alternative terminology, the lattice 
has "U(l) magnetic disorder") if these distributions are flat for Polyakov lines and large 
Wilson loops. In other words, we have U(l) symmetry iff, for any g,g' G U(l), it is true 
that 

V T [g]-V T [g'g] = (1.6) 

for Polyakov lines, and 

V c [g] - Vdg'g] ~ exp[-aArea(C)] (1.7) 

holds asymptotically for large Wilson loops. In the case that these relations are not true 
in general, but the restricted forms 

V T [g}-V T [zg} = (1.8) 

and 

Vc[g) ~ V c [zg] ~ exp[-(xArea(C)] (1.9) 

hold for any g G U(l) and z = ±1 G Z2, then we will say that the lattice has only Z2 
global symmetry (or Z 2 magnetic disorder) .0 

Generalizing to g 6 SU(N) and z 6 Zjy, for gauge-invariant loops on the unprojected lattice, eq. ( |l.8| ) 
follows from the well-known global Zn symmetry of the confined phase. 
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Inserting (|1.5| ) into ( 



Q| ) and (|1.4|), we have 



V c \e ie \ = ^ [l + 2^W q (C) cos(q6) 



V T [e W ] = ^[l + 2j2P,cos(qe)) (1.10) 



From this, we can immediately see that there is U(l) magnetic disorder iff W q (C) has 
an asymptotic area law falloff, and P q = 0, for all integer q ^ 0. On the other hand, 
if these conditions hold only for q =odd integer, then we have Z 2 magnetic disorder. It 
is an unambiguous prediction of the vortex theory that the lattice has only Z 2 magnetic 
disorder, even after abelian projection. 

The magnetic disorder induced by a monopole Coulomb gas is expected to be rather 
different from the disorder induced by vortices. A Coulombic magnetic field distribution 
will in general affect loops of any q, with q > 1 loops responding even more strongly than 
q — 1 to any variation of magnetic flux through the loop. The usual statement is that 
all integer abelian charges are confined. This statement is confirmed explicitly in D = 3 
dimensions, where it is found in a semiclassical calculation that the monopole Coulomb 
gas derived from QED3 confines all charges, with string tensions a q directly proportional 
to the charge 

a q = q<7 1 (1.11) 

This relation is also consistent with recent numerical simulations ||. 

In D=4 dimensions, an analytic treatment of monopole currents interacting via a two- 
point long-range Coulomb propagator, plus possible contact terms, is rather difficult. Nev- 
ertheless, QED4 (particularly in the Villain formulation) can be viewed as a theory of 
monopole loops and photons, and in the confined (=strong-coupling) phase there is a 
string tension for all charges q, and all P q = 0, as can be readily verified from the strong- 
coupling expansion. Confinement of all charges q is also found in a simple model of the 
monopole Coulomb gas, due to Hart and Teper B. Finally, all multiples q of electric 
charge are confined in the dual abelian Higgs model (a theory of dual superconductivity), 



and this model is known (10| to be equivalent, in certain limits, to an effective monopole 
action with long-range two-point Coulombic interactions between monopole currents. 

In the case of D = 4 abelian-projected Yang-Mills theory, we do not really know if the 
distribution of monopole loops identified on the projected lattice is typical of a monopole 
Coulomb gas. What can be tested, however, is whether or not the field associated with 
these monopoles is Coulombic (as opposed, e.g., to collimated). This is done by comparing 
observables measured on abelian-projected lattices with those obtained numerically via a 



"monopole dominance" (MD) approximation, first introduced in ref. |T2|]. The MD approxi- 
mation involves two steps. First, the location of monopole currents on an abelian-projected 
lattice is identified using the standard DeGrand-Toussaint criterion fT3"| . Secondly, a lat- 
tice configuration is reconstructed by assuming that each link is affected by the monopole 
currents via a lattice Coulomb propagator. Thus, the lattice Monte Carlo and abelian 
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projection supply a certain distribution of monopoles, and we can study the consequences 
of assigning a Coulombic field distribution to the monopole charges. 

We should pause here, to explain how the non- confinement of color charges in the adjoint 
representation is accounted for in the monopole gas or dual-superconductor pictures, in 
which all abelian charges are confined. The screening of adjoint (or, in general, j =integer) 
representations in Yang-Mills theory is a consequence of having only Z 2 global symmetry 
on a finite lattice in the confined phase. In a monopole Coulomb gas, on the other hand, 
the confinement of all abelian charges would imply a U(l) global symmetry of the projected 
lattice at finite temperature. Nevertheless, the Z 2 symmetry of the full lattice and U(l) 
symmetry of the projected lattice are not necessarily inconsistent. Let Pj M represent the 



Polyakov line in SU(2) gauge theory in group representation j. Generalizing eq. ( L4|) to 
g G SU(2), we have 

r T [g}= E Pj M xAg\ (1.12) 

3=0,1,1... 



and the fact that eq. rather than eq. ( |1.6|) , is satisfied follows from the identity 

Xj(zg) = Xj(g) for 3 = integer (1.13) 



and from the fact that 



p ym 1=0 j = half-integer . , 

j \ ^0 j = integer > 1 ' 1 



in the confined phase, due to confinement of charges in half-integer, and color-screening of 
charges in integer, SU(2) representations. 

Now, on an abelian projected lattice, the expectation value of a Yang-Mills Polyakov 
line in representation j becomes 

pJM _ £ (L15) 



where P q is the abelian Polyakov line defined in eq. (|Q). If the abelian-projected lattice 
is U(l) symmetric, then P 2m = for all m ^ 0, while Pq — 1. This means that after 
abelian-projection we have 

p YM f j = half-integer ^ lg , 

3 1 1 j — integer 

and in this way the non-zero values of j =integer Polyakov lines are accounted for, even 
assuming that the projected lattice has a global U(l)-invariance. Similarly, adjoint Wilson 
loops will not have an area law falloff on the projected lattice, as expected from color- 
screening. This explanation of the adjoint perimeter law in the abelian projection theory 
has other difficulties, associated with Casimir scaling of the string tension at intermediate 
distances (cf. ref. but at least the asympotic behavior of Wilson loops in various 



representations is consistent with U(l) symmetry on the abelian projected lattice. 
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It is just this assumed U(l) symmetry of the abelian-projected lattice, and the monopole 
Coulomb gas picture which is associated with it, that we question here. The U(l) vs. Z 2 
global symmetry issue can be settled by calculating abelian Wilson loops and/or Polyakov 
lines for abelian charges q > 1. The validity of the Coulomb gas picture can also be 
probed by calculating W q (C) and P q with and without the monopole dominance (MD) 
approximation, and comparing the two sets of quantities. 

The first calculation of q > 1 Wilson loops, in the MD approximation, was reported 
recently by Hart and Teper |J; their calculation confirms the Coulomb gas relation o q oc q 
found previously for compact QED3. In ref. ||, this result is interpreted as favoring the 
monopole Coulomb gas picture over the vortex theory. From our previous remarks, it may 
already be clear why we do not accept this interpretation. At issue is whether monopole 
magnetic fields spread out as implied by the Coulomb propagator, or whether they are 
collimated in units of ±7r. This issue cannot be resolved by the MD approximation, which 
imposes a Coulombic field distribution from the beginning. The MD approximation does, 
however, tell us that if the monopoles have a Coulombic field distribution, then the q = 2 
Wilson loop has an area law falloff, at least up to the maximum charge separation studied 
in ref. ||. The crucial question is whether the q = 2 loops computed directly on abelian 
projected lattices also have an asymptotic area law falloff, or instead go over to perimeter- 
law behavior (usually called "string-breaking") as predicted by the vortex theory. 

Here it is important to have some rough idea of where the q = 2 string is expected to 
break, according to the vortex picture, otherwise a null result can never be decisive. A 
q = 2 Wilson loop will go over to perimeter behavior when the size of the loop is comparable 
to the thickness of the vortex. It seems reasonable to assume that the thickness of a vortex 
on the abelian-projected lattice is comparable to the thickness of a center vortex on the 
unprojected lattice. From Fig. 1 of ref. ||, this thickness appears to be roughly one fermi, 
which is also about the distance where an adjoint representation string should break in 
D = 4, according to an estimate due to Michael ||16fl . The finite thickness of the vortex is 
an important feature of the vortex theory, as it allows us to account for the approximate 



Casimir scaling of string tensions at intermediate distances (cf. refs. [|T7j , [Tq1 ) . But it also 
means that at, e.g., (3 = 2.5, we should look for string breaking at around R = 12 lattice 
spacings. Noise reduction techniques, such as the "thick-link" approach, then become 
essential. 

The validity of the thick-link approach, however, is tied to the existence of a transfer 
matrix. Since the method uses R x T loops with R 3> T, one has to show that the 
potential extracted is mainly sensitive to the large separation R, rather than the smaller 
separation T, and here positivity plays a crucial rule. Since there is no transfer matrix 
in maximal abelian gauge, the validity of the thick link approach is questionable (and 
the issue of positivity is much more than a quibble, as we will see below). Moreover, 
even when a transfer matrix exists, string-breaking is not easy to observe by this method, 
and requires more than just the calculation of rectangular loops. The breaking of the 
adjoint-representation string has not been seen using rectangular loops alone, and only 
quite recently has this breaking been observed, in 2+1 dimensions, by taking account of 
mixings between string and gluelump operators |TJ], ITSf . The analogous calculation, for 
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operators defined in maximal abelian gauge, would presumably involve mixings between 
the q = 2 string and "charge-lumps"; the latter being bound states of the static abelian 
charge and the off-diagonal (double abelian-charged) gluons. 

There are, in fact, existing calculations of the q = 2 potential, by Poulis [19| and by Bali 



et al. ||20|| , using the thick-link method. String breaking was not observed, but neither did 
these calculations make use of operator-mixing techniques, which seem to be necessary for 
this purpose. In any case, in view of the absence of a transfer matrix in maximal abelian 
gauge, we do not regard these calculations as decisive. 

Given our reservations concerning the thick-link approach, we will opt in this article 
for a far simpler probe of global symmetry/magnetic disorder on the projected lattice, 
namely, the double-charged abelian Polyakov lines P<i- Any abelian magnetic vortex can 
be regarded, away from the region of non- vanishing vortex field strength, as a discontinuous 
gauge transformation, and it is this discontinuity which affects Wilson loops and Polyakov 
lines far from the region of finite vortex field strength. If the vortex flux is ±7r, the 
discontinuity will not affect even-integer g-charged Polyakov lines, and these should have a 
finite expectation value. The abelian-projected lattice then has only Z 2 global symmetry 
in the confined phase. In contrast, a monopole Coulomb gas is expected to confine all 
q charges, as in compact QED S , and the Z 2 subgroup should play no special role. In 
that case, P q = for all q. Thus, if we find that the q = 2 Polyakov line vanishes, this 
is evidence against vortex structure and flux collimation, and in favor of the monopole 
picture. Conversely, if q — 2 Polyakov lines do not vanish, the opposite conclusion applies, 
and the vortex theory is favored. 



2 Polyakov lines 

After fixing to maximal abelian gauge in SU(2) lattice gauge theory, abelian link variables 

Ufa) = diag[e iefl{:!;) , e~ ie ^ x) } (2.1) 

are extracted by setting the off-diagonal elements of link variables to zero, and rescaling 
to restore unitarity. A g-charge Polyakov line P q {x) is defined as 



P q (x) = Y[ exp[zg6 , 4 (x + nA)] 



(2.2) 



ra=l 



where is the number of lattice spacings in the time (space) directions. We can 

consider both the expectation value of the lattice average 



Pn 



and the expectation value of the absolute value of the lattice average 



jabs 



[2.3) 



(2.4) 
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Polyakov lines can vanish for q = 1, even in the deconfined phase, just by averaging over 
Z 2 -degenerate vacua, which motivates the absolute value prescription. In the confined 
phase, one then has Pi bs ~ N s 3 ^ 2 . We will find that this prescription is unnecessary for 
q = 2, and we will compute these Polyakov lines without taking the absolute value of the 
lattice average. 

For purposes of comparison, and as a probe of the monopole Coulomb gas picture, we 
also compute "monopole" Polyakov lines Pm 9 following an MD approach used by Suzuki et 
al. in ref. JHJ]. Their procedure is to decompose the abelian plaquette variable (<9 M denotes 
the forward lattice difference) 

W^W-W (2-5) 

into two terms 

f^u(x) = f^(x) + lim^x) (2.6) 

where n Mi , is an integer-valued Dirac-string variable, and —n < f < tt. One can then 
invert (|2.5|) to solve for # 4 in terms of the "photon" field-strength / , the Dirac-string 
variables n^ u , and an irrelevant U(l) gauge- dependent term. If we assume that the photon 
and Dirac-string variables are completely uncorrelated, then the Dirac-string contribution 
is given by 

e^(x) = -J2D(x,x')dln u4 (x') (2.7) 

x< 

Here D(x,x') is the lattice Coulomb propagator, and the partial derivative denotes a 
backward difference. The monopole dominance approximation is to replace 84 by O^f 
in eq. ( p. 2D , the idea being that this procedure isolates the contribution of the monopole 
fields to the Polyakov lines. The correlations between the photon, monopole, and abelian 
lattice fields will be discussed in more detail in section 3, and in an Appendix. 

2.1 Z 2 Magnetic Disorder 

In Fig. |S| we display P^ bs and P^l f° r the q — 1 lines, on a 12 3 x 3 lattice. There are no 
surprises here; we see that for = 3 there is a deconfinement transition around (3 = 2.15. 

The situation changes dramatically when we consider q = 2 Polyakov lines. Fig. |] is a 
plot of the values of P2 and Pm2, without any absolute value prescription, on the 12 3 x 3 
lattice. To make the point clear, we focus on the data in the confined phase, in Fig. 0. It 
can be seen that P2 is non-vanishing and negative in the confined phase; the data is clearly 
not consistent with a vanishing expectation value. In the MD approximation, Pmi m &y 
also be slightly negative, but its value is at least an order of magnitude smaller than P2. 
This seems to be a very strong breakdown of monopole dominance, in the form proposed 



in ref. 



In Figures |3l-|T0l we plot the corresponding data found on a 16 3 x 4 lattice. There 
is a deconfinement transition close to (3 — 2.3, and again there is a clear disagreement 
between P 2 and Pm2, with the former having a substantial no n- vanishing expectation 
value throughout the confined phase. 
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q=1 Polyakov Lines at T=3 

1 | , , , , , , r- 




2 2.05 2.1 2.15 2.2 2.25 2.3 2.35 2.4 2.45 2.5 
Beta 



Figure 5: Single charge (q = 1) Polyakov lines at time extention T=3, on the abelian 
projected lattice and in the monopole dominance (MD) approximation. 



q=2 Polyakov lines at T=3 




2 2.05 2.1 2.15 2.2 2.25 2.3 2.35 2.4 2.45 2.5 
Beta 

Figure 6: Doubly-charged (q = 2) Polyakov lines at time extension T=3, on the abelian 
projected lattice and in the monopole dominance (MD) approximation. 

The numerical evidence, for both T = 3 and T = 4, clearly favors having Z 2 , rather 
than U(l), global symmetry/magnetic disorder on the abelian projected lattice. 

2.2 Spacelike Maximal Abelian Gauge 

It is also significant that P 2 is negative. This implies a lack of reflection positivity in the 
Lagrangian obtained after maximal abelian gauge fixing, and must be tied to the fact that 
maximal abelian gauge is not a physical gauge. This diagnosis also suggests a possible 
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q=2 Polyakov lines at T=3 
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Figure 7: Magnified view of q = 2 abelian and monopole-dominance Polyakov lines at time 
extension T = 3, in the confined phase. 

q=1 Polyakov Lines at T=4 
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Figure 8: q = 1 Polyakov lines (as in Fig. [|), for time extension T=4 lattice spacings. 



cure: Instead of fixing to the standard maximal abelian gauge, which maximizes 

R = J2 y Ln°Mx)a 3 Ul(x)] 



(2.8) 



we could try to use a "spacelike" maximal abelian gauge [22[], maximizing the quantity 



E = ^^Trh[/ t (^ 3 (/|(i)] 

x k=l 



(2.9) 



which involves only links in spatial directions. This is a physical gauge. What happens 
in this case is that one disease, the loss of reflection positivity, it replaced by another, 
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q=2 Polyakov lines at T=4 
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Figure 9: q = 2 Polyakov lines (as in Fig. ||), for time extension T=4 lattice spacings. 



q=2 Polyakov lines at T=4 
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Figure 10: q = 2 Polyakov lines in the confined phase (as in Fig. [?]), for time extension 
T=4 lattice spacings. 



namely, the breaking of 90° rotation symmetry. This is illustrated in Fig. [XX], where we 
plot spacelike and timelike Polyakov lines on a 4 4 lattice, in the spacelike maximal abelian 
gauge defined above. We find that the values for double-charged Polyakov lines running in 
the time direction are much reduced in the spacelike gauge, and in fact the results shown 
appear consistent with zero. Spacelike q = 2 Polyakov lines, however, which run along 
the 1,2, or 3 lattice directions, remain negative, and in fact are larger in magnitude than 
Polyakov lines of the same length, and the same coupling, computed in the usual maximal 
abelian gauge. One therefore finds on a hypercubic lattice that 90° rotation symmetry is 
broken. 
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The spacelike Polyakov line operator creates a line of electric flux through the periodic 
lattice. The non- vanishing overlap of this state with the vacuum has, in the spacelike gauge, 
a direct physical interpretation: Since the q = 2 electric flux line cannot, for topological 
reasons, shrink to zero, a finite overlap with the vacuum means that the q = 2 flux tube 
breaks. This is presumably due to screening by double-charged (off-diagonal) gluon fields. 
The implication is that in a physical gauge, where Wilson loops can be translated into 
statements about potential energies, q = even abelian charges are not confined. 
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Spacelike Maximal Abelian Gauge, 4**4 Lattice 
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Figure 11: Breaking of 90° rotation symmetry in spacelike maximal abelian gauge, as seen 
from comparing timelike Polyakov lines (lattice 4 direction) with spacelike Polyakov lines 
(lattice 1,2, or 3 directions) on a 4 4 lattice. 

The fact that q = even charges are unconfined, together with the positivity property 
in spacelike maximal abelian gauge, leads to the conclusion that timelike q = 2 Polyakov 
lines P2 are positive and non-zero, although the data points for timelike P2 shown in Fig. 
|TT| , which appear to be consistent with zero, do not yet support such a conclusion. It must 
be that the value of P 2 in this gauge is simply very small, and much better statistics are 
required to distinguish that value from zero. To get some idea of the difficulty involved, 
let us suppose that the magnitude of the timelike abelian line P2 on the projected lattice 
is comparable to the magnitude of the gauge- invariant Polyakov line P a dj, in the adjoint 
representation of SU(2), on the full, unprojected lattice. To leading order in the strong- 
coupling expansion, P a ^ on a lattice with extension T in the time direction is given by 

P «* = 4 (f) ( 2 - 10 ) 

This equals, e.g., 0.00156 for T = 2 at (3 = 1.5; quite a small signal considering that T 
is only two lattice units. The obvious remedy is to increase (3, but then one runs into a 
deconfinement transition at (3 = 1.8. We can move the transition to larger (3 by increasing 
T, but of course increasing T again causes the signal to go down. 
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The best chance to extract a signal from the noise is to choose a value of (3 which is fairly 
close to the deconfmement transition (but still in the confined phase), and to generate very 
many configurations. Here are the results obtained at T = 2 lattice spacings and [3 — 1.7, 
coming from 5000 configurations separated by 100 sweeps on a 2 x 8 3 lattice: 

P ad] = 0.00447(23) 
P 2 = 0.00241(52) (2.11) 

The result for the adjoint line is consistent with the strong-coupling prediction of P a dj = 
0.00426. The abelian line P 2 is non-zero, positive, and comparable in magnitude to P a dj, 
although it must be admitted that the errorbar is uncomfortably large. The corresponding 
values for T = 3 and (3 = 2.11, obtained from 5000 3 x 8 3 lattices, are 

p adj = 0.00338(28) 
P 2 = 0.00124(42) (2.12) 

Again P 2 is non-zero, although the errorbar is still too large for comfort. 

Clearly the evaluation of the timelike P 2 line in spacelike maximal abelian gauge is cpu- 
intensive, and our results for this quantity must be regarded as preliminary. Nevertheless, 
these preliminary results are consistent with the conclusion previously inferred from the 
spacelike lines: In a physical gauge, the q = 2 charge is screened, rather than confined, 
and we have Z 2 , rather than U(l), magnetic disorder on the abelian-projected lattice. 



3 The "Photon" Contribution 

Suppose we write the link angles O^x) of the abelian link variables as a sum of the link 
angles dff(x) in the MD approximation, plus a so-called "photon" contribution 9^ l (x), i.e. 

ef{x) = e»{x)-e™{x) (3.1) 



It was found in refs. [|T2, ZT|, 23"| that the photon field has no confinement properties at all; 
the Polyakov line constructed from links U ^ = exp[i<?2 ] is finite, and corresponding Wilson 
loops have no string tension. Since 9*f would appear to carry all the confining properties, 
a natural conclusion is that the abelian lattice is indeed a monopole Coulomb gas. 

To see where this reasoning may go astray, suppose we perversely add, rather than 
subtract, the MD angles to the abelian angles, i.e. 



M j 

0f{x) + K%(x) (3.2) 



e'^x) = e,(x) + e™(x 



in effect doubling the strength of the monopole Coulomb field. It is natural to expect a 
corresponding increase of the string tension, and of course Pi = should remain true. 
Surprisingly, this is not what happens; doubling the strength of the monopole field in fact 
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removes confinement .Q Some results for Pi are shown in Table |1| Here we have computed 
the vev of P\(x) without taking the absolute value of the lattice sum (i.e. we use eq. 
( j2.3| ) rather than eq. (IO)), and we find that Pi is finite and negative in the additive 
configurations. The additive configuration 9' is far from pure-gauge, and the vev of P\ is 
correspondingly small. Nevertheless, Pi is non-zero, so adding the monopole field in this 
case actually removes confinement. Clearly, the interplay between the MD and "photon" 
contributions is a little more subtle than previously supposed. 



T 


(3 


Pi line 


3 


1.8 


-0.0299 (20) 


3 


2.1 


-0.0405 (10) 


4 


2.1 


-0.0134 (10) 



Table 1: Single-charged abelian Polyakov lines, computed in the additive 9' = 9 + 9 M 
configurations, in which the strength of the MD contribution is effectively doubled. 



To understand what is going on, we return to the concept of the holonomy probability 
distribution 

N T 

5 e ie , J] exp[i9 A (x + nA)] ) (3.3) 

n=l 



V(6) = V T [e 
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Pr[e i61 ] is the probability density for the U(l) group elements on the group manifold. 
However, since the group measure on the U(l) manifold is trivial (i.e. dQ), it is not hard to 
see that V(9)d9 is interpreted as the probability that the phase of an abelian Polyakov line 
lies in the interval [9, 9 + 69]. In a similar way, replacing 9± by 9^ or 9% h on the rhs of the 
above equation, we can define the probability distributions V(9 AI ) and V(9 ph ), respectively, 
for the phases of monopole and photon Polyakov lines. All of these distributions have 2n- 
periodicity, and are invariant under 9 — > —9, reflections, so we need only consider their 
behavior in the interval [0, 7r]. Without making any further calculations, it is already 
possible to deduce something about the shape of V(9): 

• Since all P q are small, V(9) is fairly flat. 

• Assuming Z 2 symmetry, V(9) is symmetric, in the interval [0, 7r], around 9 = | . 

• Since P 2 is negative, V(9) should be larger in the neighborhood of 9 = ~ than in the 
neighborhood of 9 = or 9 = n. 

From these considerations, we deduce that V(9) looks something like Fig. [12]. Similarly, 
since PMq ~ in the MD approximation, we conclude that there is very nearly U(l) 
symmetry in this approximation, and V(9 M ) is almost flat, as in Fig. ITS. 



3 We have already noted that in the absence of a transfer matrix, the term "confinement of abelian 
charge" must be used with caution. In this section, the phrase "confinement of charge q" is just taken to 



mean ll P q = 0" . 
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Abelian Link Angle Distribution 




71/2 
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Figure 12: Conjectured V{9) distribution of the Polyakov phase angle, based on Z% global 
symmetry and P2 < 0. 

Monopole Angle Distribution 



X3 
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71/2 

Monopole Angle 6 M 



Figure 13: V(9 ) in the MD approximation, assuming perfect global U(l) symmetry. 



Now if the link angles 9 fl (x) and 9^(x) are correlated to some extent, then the differ- 
ence 8 ph (x) between these variables is not random, but has some non-uniform probability 
distribution as illustrated in Fig. [T3| Since the Fourier cosine components of V(9 ph ) are 
typically non-zero, it follows that the photon field, by itself, has no confinement property. 
The crucial point is that by subtracting 9 , the Z2 symmetry of 9 is broken due to the 
correlation between 9 and 9 M . It is interesting to note that even if V(9 M ) were neither U(l) 
nor Z 2 symmetric, i.e. if we imagine that the ^-configurations confine but the MD contri- 
butions do not, a correlation between the 9 and 9 M would still be sufficient to break the 
Z2 symmetry of the difference configuration 9 ph . As a result, subtracting the non-confining 
9 M from the confining 9 would still remove confinement. 

In the center vortex picture, vortex fields supply the confining disorder, but of course 
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Photon Angle Distribution 



o 



I 1 

-71 (J 31 

Photon Angle e ph 

Figure 14: The "photon" angle distribution V(6 ph ), resulting from a correlation between 
9and9 M . 

this does not at all exclude a correlation of the MD variables 9 M with the 9 variables. 
According to the arguments in the Introduction, monopoles lie along vortices as shown in 
Fig. 1 (further evidence is given in the next section), and this correspondence will certainly 
introduce some degree of correlation between magnetic flux on the abelian-projected lattice, 
and magnetic flux in the MD approximation. Roughly speaking, one can say that the 
confining flux has the same magnitude on the abelian and MD lattices, only it is distributed 
differently (collimated vs. Coulombic). However, according to the center vortex picture, 
there must also be some correlation between 9 p Mx) and 9jf(x); this is necessary to convert 
the long-range monopole Coulomb field into a vortex field, and to break the U(l) symmetry 
of the MD lattice down to the Z 2 symmetry of the vortex vacuum. 

In numerical simulations performed at (3 = 2.1 and T = 3 on a 3 x 12 3 lattice, we do, 
in fact, find a striking correlation between 9 ph and 9 M : The average "photon" angle 9 ph 
tends to be positive for 9 M e [0, ~], and negative for 9 M 6 [| ,ir). Computing the average 

photon angle W h in each monopole angle quarter-interval, we find 

W k _ j 0.027(4) for 9 M e [0, f ] 

~\ -0.027(4) for^ M G[f,vr] ^' 4j 
This result, combined with the results displayed in Table [IJ, raises two interesting questions: 

1. How is the correlation between 9 ph and 9 M , found above in ( |3.4| ), related to the 
remaining Z 2 global symmetry of V(9)\ and 

2. Why is Pi negative in the additive configurations of eq. (|3.2j) ? 

To shed some light on these issues, we begin by defining 9 ph (9 M ) as the average value 
of 9 ph at fixed 9 M , and then make the drastic approximation of neglecting all fluctuations 
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of 9 ph at fixed 9 M around its mean value. This amounts to approximating the vev of any 
periodic function F{9) 

< F >= ^ d9 F{9)V{9) (3.5) 

by 

<F>= T d9 M ^F(9 M + 9 P \9 M )) (3.6) 

where we have used the fact that the probability distribution for 9 M is (nearly) uniform. 
The accuracy of this approximation depends, of course, on the width of the probability 
distribution for 9 ph at fixed 9 M , and on the particular F{9) considered. Here we are only 
concerned with certain qualitative aspects of phase angle probability distributions, and 
hopefully the neglect of fluctuations of 9 ph around the mean will not severely mislead us. 

With the help of the approximation (|3.6|), we can answer the two questions posed above. 
In this section we will only outline the argument, which is presented in full in an Appendix. 

The function '9 ph {9 M ) maps the variable 9 M E [— it, tt], which has a uniform probability 
distribution in the interval, into the variable 9 G [— vr,7r], where 

e = 9 M + 9 p \9 M ) (3.7) 
The non-uniform mapping induces a non-uniform probability distribution for the variable 

1 d9 M 



V{9) 




d£^ 
d9 



(3.8) 



which we identify with V{9) in the approximation ( |3.6|) . Since V{9) is peaked at 9 = ±|, 

it follows that d9 ph ' jd9 is minimized at 9 — ±|. 

Global Z<i symmetry implies that P q = for q = odd. Then, from eq. ( |1.10| ) we have 



V(ir-9) = V{9) 
V{-ti-9) = V(9) 

V(-9) = V(9) (3.9) 

~nP h ( qM\ _ 



From these relationships, eq. (|3.8|) , and the fact (shown in the Appendix) that 9 (—9 
-9 ph (9 M ), we find that 

9 ph [9]=9 ph [9 M (9)] (3.10) 

is an odd function with respect to reflections around 9 = 0, ±|. Defining as the average 
9 ph in the quarter interval 9 M £ [0, |], and 9^ as the average 9 ph in the quarter-interval 
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9 M e [f ,7r], we have 

= £ r /2 d9 M r h (e M ) 

IT Jo 

ri 2 - - -oh - 

JO 

= 4 T d9 V(n -9)9 ph [n -9} 

Jtt/2 

= -4 r 

Jtt/2 

= _£ r ^ gph/^ 

7T Jtt/2 

= -SS (3.H) 

which explains, as a consequence of global Z 2 symmetry, the equal magnitudes and opposite 
signs found in eq. (|3.4j ). This answers the first of the two questions posed above. 

For expectation values of Polyakov phase angles in the additive configuration, we have 

< F >= r d& F{9')V'{9') (3.12) 



where V(9') is the probability distribution for the Polyakov angles of the additive config- 
uration 9' = 9 + 9 AI . Again neglecting fluctuations of 9 ph around the mean r\9 M ), and 
changing variables to 9 = 9 + 9 M we have 



< F > = r d9 M —F{29 M + r h (9 M )) 



r2n , i dQ M 

/ d 9' —^=rF(9' 

J2n 27T dff 



2vr d6 

d9 ! -^=rF(9) (3.13) 

7T d9 



where the 27r-periodicity of the integrand was used in the last step. Then the induced 
probability distribution in the 9 variable is 

vr d9 

7i y d9 J \ d9 J y ' 

As shown in the Appendix, 9 = ±^ corresponds to 9 = ±7r, and the assumed single- 
valuedness of 9 M (9) requires d9 ph /d9 < 1. In that case, since d9 ph /d9 is minimized at 
9 = ±~, it follows that V(9 ) has a peak at 9 = ±7r. Given that the 9' distribution is 
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peaked at ±7r, as in Fig. the n = 1 coefficient in the cosine series expansion of this 
distribution (which by definition is Pi) is evidently negative, answering the second of the 
two questions posed below ( |3.4| ). Confinement is lost because the Z2 symmetry of the 9 
distribution has been broken. 

Additive Config. Angle Distribution 



■8 



71/2 

Additive Angle 6' 



Figure 15: Polyakov phase distribution V'{9') in the additive configuration 6' = 6 + . 



The correlation that exists between the monopole and photon contributions in abelian- 
projected SU(2) gauge theory implies that these contributions actually do not factorize in 
Polyakov lines and Wilson loops, in contrast to the factorization which occurs in compact 
QED in the Villain formulation. In fact, the terminology "photon contribution" used to 
describe 9 ph is really a little misleading. The field O p J l {x) is best described as simply the 



difference O^x) 







M , 



X 



between the abelian angle field and the MD angle field. It is not 



correct to view 6 ph as a purely perturbative contribution, since the correlation that exists 
between 6 ph and Off , which breaks U(l) down to an exact Z2 remnant symmetry, clearly 
has a non-perturbative origin. 



Finally, in Figs. |16l|18t we show some histograms for the actual probability distributions 
V(6),V(6 M ),V(6 ph ), obtained on a 3 x 12 3 lattice at (3 = 2.1. V(0) and V(6 M ) are shown 
on the [0, tt] half-interval, while V(6 ph ) is displayed on the full [— tt, tt] interval. The height 
of the histogram is the probability for \0\, \6 M \, 6 ph to fall in each interval. It is clear that 



these numerical results agree with the conjectured behavior in Figs. 12-14. 



4 Field Collimation 

Although the finite VEV of q = 2 Polyakov lines is a crucial test, it is also useful to ask 
whether the collimation of confining field strength into vortex tubes can be seen more 
directly on the lattice. 

In the most naive version of the monopole Coulomb gas, the monopole field is imagined 
to be distributed symmetrically, modulo some small quantum fluctuations, around a static 
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Figure 16: Histogram of V{9) in the interval [0, 7r], for T = 3 links in the time direction at 
(3 = 2.1. Note the symmetry around 9 = tt/2, which is associated with Z 2 symmetry. 
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Figure 17: Histogram of V(6 M ), for Polyakov lines in the MD approximation. Time extent 
is T = 3 at = 2.1. 



monopole. In this section we will find that the field around the position of an abelian 
monopole, as probed by SU(2)-invariant Wilson loops, is in fact highly asymmetric, and 
is very strongly correlated with the direction of the center vortex passing through the 
monopole position. Some of these results, for unit cubes around monopoles, have been 
reported previously in ref. H, but are included here for completeness. The results for 3- 
and 4-cubes around static monopoles are new. 

To circumvent the Gribov copy issue, we work in the "indirect" maximal center gauge 
introduced in ref. ||24j| , and locate monopole and vortex positions by projections (abelian 
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Figure 18: Histogram of V(8 ph ), for the "photon" Polyakov lines, on the full interval 
[— 7T, 71"]. Time extent is T — 3 at /3 — 2.1. 



and center, respectively) of the same gauge-fixed configuration. Indirect maximal center 
gauge is a partial fixing of the U(l) gauge symmetry remaining in maximal abelian gauge, 
so as to maximize the squared trace of abelian links. The residual gauge symmetry is 
Z2- The excitations of the center projected lattice are termed "P- vortices," and have been 
found to lie near the middle of thick center vortices on the unprojected lattice (cf. H). 



4.1 Monopole-Antimonopole Alternation 

According to the argument depicted in Figs. [|-[|, at any fixed time the monopoles found 
in abelian projection should lie along vortex lines, with monopoles alternating with an- 
timonopoles along the line. To test this argument, we consider static monopoles (associ- 
ated with timelike monopole currents) on each constant time volume of the lattice. Each 
monopole is associated with a net ±27r magnetic flux through a unit cube. In numerical 
simulations performed at (3 = 2.4, we find that almost every cube, associated with a static 
monopole, is pierced by a single P-vortex line. Only very small fractions are either not 



pierced at all, or are pierced by more than one line, with percentages shown in Fig. [19. 

P- vortices are line- like objects on any given time slice of the lattice.^ About 61% of these 
vortex lines have no monopoles at all on them. We find that 31% contain a monopole- 
antimonopole pair. The remaining 8% of closed vortex lines have an even number of 
monopoles + antimonopoles, with monopoles alternating with antimonopoles as one traces 
a path along the loop. This is exactly the situation sketched in Fig. [11. Exceptions to 
the monopole- antimonopole alternation rule were found in only 1.2% of loops containing 
monopoles. In every exceptional case, a monopole or antimonopole was found within one 



4 It should be noted that vortices are surface-like objects in D=4 dimensions, so different closed loops 
on a given time slice may belong to the same P-vortex surface. 
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lattice unit of the P-vortex line which, if counted as lying along the vortex line, would 
restore the alternation. 




3 % 93 % 4 % 

No vortex 1 vortex >1 vortex 

Figure 19: Percentages of monopole cubes pierced by zero, one, or more than one P-vortex 
lines, at /3 = 2.4. 



4.2 Field Collimation on 1-cubes 

We define vortex limited Wilson loops W n (C) as the expectation value of Wilson loops 
on the full, unprojected lattice, subject to the the constraint that, on the projected lat- 
tice, exactly n P- vortices pierce the minimal area of the loop (cf. 0). We employ these 
gauge-invariant loop observables to probe the (a)symmetry of the color field around static 
monopoles, again at (3 = 2.4. 

Consider first, on a fixed time hypersurface, the set of all unit cubes which contain 
one static monopole, inside a cube pierced by a single P-vortex line. This means that two 
plaquettes on the cube are pierced by the vortex line, and four are not. The difference S 
between the average plaquette Sq on the lattice, and the plaquette on pierced/unpierced 
plaquettes of the monopole cube 

S = S -< IttPUU^)^ face > (4.1) 



is shown in Fig. EO. For comparison, we have computed the same quantities in unit cubes, 



pierced by vortices, which do not contain any monopole current. 

It is obvious that the excess plaquette action associated with a monopole is extremely 
asymmetric, and almost all of it is concentrated in the P-vortex direction. Moreover, the 
action distribution around a monopole cube is not very different from the distribution on 
a cube pierced by a vortex, with no monopole at all inside. The two distributions are 
even more similar, if we make the additional restriction to "isolated" static monopoles; i.e. 
monopoles with no nearest-neighbor monopole currents. The excess action distribution for 
isolated monopoles, again compared to zero-monopole one-vortex cubes, is shown in Fig. 
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Figure 20: Excess plaquette action distribution on a monopole cube pierced by a single 
P-vortex. For comparison, the excess action distribution is also shown for a no-monopole 
cube pierced by a P-vortex. 
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Figure 21: Same as Fig. 20 for isolated monopoles 



4.3 Field Collimation on Larger Cubes 

Finally we consider cubes which are N lattice spacings on a side, in fixed-time hypersur- 
faces, having two faces pierced by a single P-vortex line and the other four faces unpierced. 
Again we restrict our attention to cubes containing either a single static monopole, or no 
monopole current. Each side of the cube is bounded by an N x N loop. Let 



Wf (N, N) 
W™(N,N) 
W?(N,N) 



1-vortex loops, bounding a monopole N-cube 

0- vortex loops, bounding a monopole N-cube 

1- vortex loops, bounding a 0-monopole N-cube 
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Wq(N, N) = 0- vortex loops, bounding a 0-monopole N-cube 



(4.2) 



denote the expectation value of N x N Wilson loops on 0/1-vortex faces of 0-monopole/ 1- 
monopole N-cubes. As a probe of the distribution of gauge- invariant flux around an N-cube, 
we compute the fractional deviation of these loops from Wq(N, N) (which has the largest 
value) by 



AM 

4° 



m W§(N,N)-W$(N,N) 



W°{N,N) 
W§(N,N)-W§ tl (N,N) 



and of course Aq = by definition. 



Wg(N,N) 



(4.3) 
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Figure 22: Gauge field collimation on a 3-cube pierced by a single P-vortex, for both one- 
monopole and no-monopole 3-cubes. A is the fractional deviation of W(3, 3) on a cube 
face from the no-monopole, no- vortex value (3, 3). 



The results, for A^ = 3 cubes and N = A cubes, are displayed in Figs. |22"| - |2~3"| , with 
the actual values for the various loop types listed in Table As with the excess-action 
distribution around a 1-cube, shown in Figs. pO|-|2T|, it is clear that gauge-invariant Wilson 
loop values are distributed very asymmetrically around a cube, and are strongly correlated 
with the direction of the P-vortex line. The presence or absence of a monopole inside the 
Af-cube appears to have only a rather weak effect on the value of the loops around each 
side of the cube; the main variation is clearly due to the presence or absence of a vortex 
line piercing the side. Obviously, the strong correlation of loop values with vortex lines, 
and the relatively weak correlation of loop values with monopole current, fits the general 
picture discussed in the introduction, of confining flux collimated into tubelike structures. 
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Figure 23: Same as Fig. 22, for 4-cubes. 



(N)-cube 


Wf (N, N) 


W^(N,N) 


W?(N,N) 


W»(N,N) 


1- cube 

2- cube 

3- cube 

4- cube 


0.3349 (10) 
0.0599 (8) 
0.0045 (7) 

-0.0030 (7) 


0.6007 (5) 
0.2426 (5) 
0.0825 (5) 
0.0255 (7) 


0.4340 (5) 
0.1101 (6) 
0.0179 (5) 
0.0016 (6) 


0.6472 (2) 
0.2587 (3) 
0.0864 (4) 
0.0268 (5) 



Table 2: Wilson loops on the faces of one-monopole (Wjff) and zero monopole (W^) 
N-cubes. The 0, 1 subscript refers to the number of P-vortices piercing the loop. 



It seems clear that the gauge field-strength distribution around monopoles identified in 
maximal abelian gauge is highly asymmetric, and closely correlated to P- vortex lines, as one 
would naively guess from the picture shown in Fig. 1. This fact, by itself, doesn't prove the 
vortex theory. One must admit the possibility that if the monopole positions were somehow 
held fixed, fluctuations in the direction of vortex lines passing through those monopoles 
might restore (on average) a Coulombic distribution. We therefore regard the findings in 
this section as fulfilling a necessary, rather than a sufficient, condition for the vortex theory 
to be correct. According to the vortex theory, confining fields are collimated along vortex 
lines, and this collimation should be visible on the lattice in some way, e.g. showing up 
in spatially asymmetric distributions of Wilson loop values. This strong asymmetry is, in 
fact, what we find on the lattice. 
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5 The Seiberg-Witten confinement scenario and 
Z2-fluxes 



The beautiful description of confinement by monopoles when N=2 supersymmetric Yang- 
Mills theory is softly broken to N— 1 supersymmetric Yang-Mills theory has enhanced 
the impression that the mechanism for confinement in pure Yang-Mills theory is monopole 
condensation. Above we have provided evidence that it is in fact Z2 fluctuations which 
are responsible for confinement in Yang-Mills theory even on abelian-projected lattices; 
see refs. || |j] for the evidence on the full, unprojected lattices. As regards the N = 2 
supersymmetric Yang-Mills theory softly broken to iV= 1, we would simply argue that 
the approximate low-energy effective theory of Seiberg and Witten is not able describe all 
aspects of confinement at sufficiently large distance scales. In particular, the low-energy 
effective theory cannot explain the perimeter law of double-charged Wilson loops. 

It may be useful here to make a distinction between the full effective action, obtained 
by integrating out all massive fields, and the "low- energy" effective action, which neglects 
all non-local or higher- derivative terms in the full effective action. The Seiberg-Witten 
calculation was aimed at determining the low-energy effective action of the softly broken 
N=2 theory. However, in a confining theory, non-local terms induced by massive fields 
can have important effects at long distances. 

The relation between the Seiberg-Witten theory and pure Yang-Mills theory in four 
dimensions has many similarities to the relation between the Georgi-Glashow model and 
pure Yang-Mills theory in three dimensions.^ In both theories the presence of a Higgs field is 
of utmost importance for the existence of monopoles. In the Georgi Glashow model it is the 
existence of a monopole condensate which is responsible for the mass of the (dual) photon 
and for confinement of the smallest unit of electric U(l) charge. However, as emphasized 
in ||, the effective low energy Coulomb gas picture of monopoles does not explain the 
fact that double-charged Wilson loops follow a perimeter law rather than an area law. 
The reason is that in a confining theory there is not equivalence between low- energy and 
long-distance physics. At sufficiently long distance it will always be energetically favorable 
to excite charged massive W^-fields and screen q = even external charges, thus preventing 
a genuine string tension between such charges. At these large scales, a description in 
terms of U(l)-disordering configurations (the monopole Coulomb gas) breaks down, and a 
description in terms of Z 2 disorder must take over. As shown in || the range of validity of 
the monopole Coulomb gas picture decreases as the mass of the H^-field decreases, and, in 
the limit of an unbroken SU(2) symmetry, confinement can only be described adequately 
in terms of Z2 fluctuations. 

In the Seiberg-Witten theory, N= 2 supersymmetry ensures that the Higgs vacuum is 
parameterized by an order parameter u = (Tr0 2 ) corresponding to the breaking of SU(2) 
to U{\). For large values of u we have a standard scenario: at energy scales fi ^> -Ju all 

5 This analogy has been used to explain the measured (pseudo)-critical exponents in four-dimensional 
lattice f/(l)-theory from the properties of N—2 supersymmetric Yang-Mills theory and its possible sym- 
metry breakings to N— 1 and N— (26) . 
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field theoretical degrees of freedom contribute to the /3-function, which corresponds to the 
asymptotically free theory. For energies lower than y/u only the U(l) part of the theory is 
effective. In these considerations the dynamical confinement scale 

A^ 2 = Axp(-87r7<?( M ) 2 ) 

obtained by the one-loop perturbative calculation plays no role. The remarkable observa- 
tion by Seiberg and Witten was that even when u < A^ =2 , where one would naively expect 
that non-Abelian dynamics was important, the system remains in the U(l) Coulomb phase 
due to supersymmetric cancellations of non-Abelian quantum fluctuations. As u decreases 
the effective electric charge associated with the unbroken U(l) part of SU(2) increases, 
while the masses of the solitonic excitations which are present in the theory, will decrease. 
Dictated by monodromy properties of the so-called prepotential of the effective low en- 
ergy Lagrangian, the monopoles become massless at a point u ~ A^ =2 where the effective 
electric charge has an infrared Landau pole and diverges. However, in the neighborhood 
of u ~ A^ =2 , this strongly coupled theory has an effective Lagrangian description as a 
weakly coupled theory when expressed in terms of dual variables, namely a monopole 
hyper-multiplet and a dual photon vector multiplet. The perturbative coupling constant 
is now go = Air/g and the point where monopoles condense corresponds to go — 0. 

A remarkable observation of Seiberg and Witten is that the breaking of N= 2 to N= 1 
supersymmetry by adding a mass term superpotential will generate a mass gap, originating 
from a condensation of the monopoles. By the dual Meissner effect this theory confines the 
electric U(l) charge at distances larger than the inverse N=2 symmetry breaking scale. In 
terms of the underlying microscopic theory it is believed that the reduction of symmetry 
from N= 2 to N= 1 allows excitations closer to generic non-supersymmetric "confinement 
excitations" , but that the soft breaking ensures that the theory is still close enough to the 
N=2 to remain an effective U(l) theory. Thus we see that in the Seiberg- Witten scenario 
we can, by introducing the mass term superpotential, describe a U(l) confining-deconfining 
transition from the N= 2 Coulomb phase to the N= 1 confining phasef]. 

But precisely as for the Georgi-Glashow model, the monopole condensate picture for 
the N — 1 confining theory is incomplete in the sense that it cannot describe the obvious 
fact that double-charged Wilson loops will have a perimeter law rather than an area law. 
Clearly, q = even external charges can be screened by the massive charged W-fields in 
the softly broken N = 2 supersymmetric Yang-Mills theory, a fact which has profound 
implications for the large-scale structure of confining fluctuations. But neither the non- 
local effects of the W^-fields, nor the VT-fields themselves, appear in the low energy effective 
action, a fact which illustrates once more that long distance physics is not captured by the 
(local) low energy effective action in a confining theory. 

6 The breaking down to N = has been analyzed in a number of papers [^7|,|2£|, where soft breaking 
via spurion fields of N= 1 and N=2 supersymmetric gauge theories are discussed (see also |29(| ). These 
models are somewhat closer to realistic models for QCD confinement, but the conclusions are, from our 
perspective, the same as for the original Seiberg- Witten model, so we will not discuss these models any 
further. 
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6 Discussion 



A point which was stressed both in ref. || and in the last section (see also [p5H), and 
which is surely relevant to the results reported here, is that charged fields in a confining 
theory can have a profound effect on the far-infrared structure of the theory, even if those 
fields are very massive. As an obvious example, consider integrating out the quark fields 
in QCD, to obtain an effective pure gauge theory. This effective pure gauge theory does 
not produce an asymptotic area law falloff for Wilson loops, which means that confining 
field configurations are somehow suppressed at large distances. A second example is the 
Georgi-Glashow model in D=3 dimensions (GG3), as discussed in ref. [|]. In this case 
the W-bosons are massive, and if their effects at large scales are simply ignored, then 
the model would be essentially equivalent to the theory of photons and monopoles, i.e. 
a monopole Coulomb gas, analyzed many years ago by Polyakov M. In the monopole 
Coulomb gas, all multiples of the elementary electric charge are confined; but this is not 
what actually happens in the Georgi-Glashow model. The reason is that W-bosons are 
capable of screening even multiples of electric charge, which means that even-charge Wilson 
loops fall only with a perimeter law, and even-charge Polyakov lines have finite vacuum 
expectation values in the confined phase. If we again imagine integrating out the W and 
Higgs fields, then the effective abelian theory confines only odd multiples of charge, the 
global symmetry is Z 2 , rather than U(l), and the theory is clearly not equivalent to either a 
monopole Coulomb gas, or to compact QED 3 . If one asks: how can the effective long-range 
theory, which involves only the photon field, be anything different from compact QED 3 , 
the answer is that the integration over W and Higgs fields produces non-local terms in the 
effective action. We note, once again, that charged fields in a confining theory have very 
long-range effects. The fact that these fields are massive does not imply that they can 
only lead, in the effective abelian action, to local terms, or that the non-local terms can 
be neglected at large scales. These remarks also apply to the Seiberg-Witten model, as 
discussed in the last section. 

In this article we have concentrated largely on a third example: abelian-projected 
Yang-Mills theory in maximal abelian gauge. Calculations on the abelian-projected lattice 
can be always regarded as being performed in an effective abelian theory, obtained by 
integrating out the off-diagonal gluon fields (and ghosts) in the given gauge. It is often 
argued that the off-diagonal gluon fields are massive, and therefore do not greatly affect 
the long-range structure of the theory. The long-range structure, according to that view, is 
dominated exclusively by the diagonal gauge fields (the "photons" ) and the corresponding 
abelian monopoles, which together are equivalent to a Coulomb gas of monopoles (D=3) or 
monopole loops (D=4). Then, since only abelian fields are involved, the global symmetry 
of the effective long-range theory is expected to be U(l), and all multiples of abelian charge 
are confined. We have seen that reasoning of this sort, which neglects the long-distance 
effects of massive charged fields, can lead to erroneous conclusions. In fact, we have found 
that on the abelian projected lattice: 

• Confinement of all multiples of abelian charge does not occur on the abelian-projected 
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lattice; charge q = 2 Polyakov lines have a non-zero VEV. 

• As a result, the global symmetry of the abelian-projected lattice is at most Z2, rather 
than U(l). 

• Monopole dominance breaks down rather decisively, at least when applied to charge 
q = 2 operators. 

• The distribution of Wilson loop values is highly asymmetric on an N-cube. There is 
a very strong correlation between loop values and the P-vortex direction, but only 
a rather weak correlation with the presence or absence of a static monopole in the 
N-cube. 

In addition, in the usual maximal abelian gauge, there is a breakdown of positivity, which 
is surely due to the absence of a transfer matrix in this gauge. The loss of positivity can 
be avoided (at the cost of rotation invariance) by going to a spacelike maximal abelian 
gauge, where we again find q = 2 string-breaking and deconfinement. The picture of a 
U(l)-symmetric monopole Coulomb gas or dual superconductor, confining all multiples 
of the elementary abelian charge, is clearly not an adequate description of the abelian- 
projected theory at large distance scales. On the other hand, the results reported here 
fit quite naturally into the vortex picture, where confining magnetic flux on the projected 
lattice is collimated in units of ±7r. 

The center vortex theory has a number of well-known (and gauge- invariant) virtues. 
In particular, the vortex mechanism is the natural way to understand, in terms of vacuum 
gauge-field configurations, the screening of color charges in zero N-ality representations, 
as well the loss of global symmetry in the deconfinement phase transition [[Mj. Center 
vortex structure is visible on unprojected lattices, through the correlation of P-vortex loca- 
tion with gauge- invariant observables [Q,|3J] . The evidence we have reported here, indicating 
vortex structure on large scales even on abelian-projected lattices, increases our confidence 
that center vortices are essential to the mechanism of quark confinement. 
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A Appendix 

In this Appendix we present the detailed argument, outlined in section 3, that of 1 = —9' 
and Pi < in the additive 9' = 9 + 9 M configurations. 
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The approximation used here is to ignore, at fixed 9m, the fluctuctions of 9 ph around 

eriodic function F{9) of the Polyakov phase 9 

d9 F(9)V{9) (A.l) 



the mean valued (9 M ); i.e. the vev of any periodic function F(9) of the Polyakov phase 9 



< F >- 



is approximated by 



< F >= f d9 M —F(9 M + 9 ph (9 M )) (A.2) 

where the factor of l/2ir corresponds to the uniform probability distribution for 9 M . The 
mean value 9 ph (9 M ) is defined as 

1 



9 ph (9 M ) 



JQM 



JQM 



where 



DO^x) 8 



Pi(x) 



Pmi(x) 



DO^x) arg [Pi{x)e 
P M i(x),e ie ' 



8 



P M1 {x),e 



(A.3) 



{ [ exp[z'0 4 (f + n4)] 

n=l 

N T 

J] exp[z0f (x + n4)] 

n=l 



(A.4) 



are Polyakov lines in the abelian and MD lattices, and S e ff is the effective abelian action, 
obtained after integrating out all off-diagonal gluons and ghost fields. Due to translation 



invariance, 9 



QP h (QM\ 



does not depend on the particular spatial position x chosen in (|A.3|) . 



From its definition, ~9 ph {9 M ) is obviously periodic w.r.t. 9 AI 



9 



M 



odd function of 9 M , i.e. 



9 ph (9 M ) 



~7ph i 



-9 IVJ 



2tt. It is also an 



(A.5) 



This is derived by first noting that the 9jf(x) link angles are functions of the 9 ll (x) link 
angles according to eqs. (|2.5|) - (|2.7|) , and that 9^(x) — > —9^{x) under the transformation 
9fi(x) — > —9 fl (x). Then, making the change of variables 9 /Ji (x) — ► —9 fl (x) in the integral 
( |A.3| ), we have 



JQM 



D9^x) 8 



DO^x) 8 



Pmi(x), e 
P M i(x),e 



-8 M 



(A.6) 



and 



r h (9 M ) 



1 



/^(x)arg (P^{x)e 



Pmi( x )-> e% 



e Seff 



-e M 



D9 ll (x)(-l) x arg (P 1 {x)e 



PM^ (x),e~ 



-9 ph (-9 M ] 



(A.7) 
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The fact that W (9 M ) is an odd function of 9 M , combined with 27r-periodicity, gives us 

^(vr) = 6 ph (-7i) = 9 ph (0) = (A.8) 
We now define the variable 

6(6 M ) = e M + e ph (e M ) (A.9) 



which is the average Polyakov phase at fixed 9 M . It will be assumed that 9(9 ) is a single- 
valued function of 9 M . Eq. ( A.9Q can then be inverted to define 9 M implicitly as a function 
of 9 

e M{e) = e-e ph [e M (e)} (a.io) 

and it will be convenient to introduce the notation 



f h [e] =r h [9 Mrs 



Applying the change of variable (|A.9|) to eq. (|A.2|) , we have 

< F >= 



<W -^Cf(9) 



2tt d9 



(All) 



(A.12) 



where, from eqs. ( |A.8|) and (|A.9| ), we see that the limits of integrations are unchanged. 
Comparing ( |A.12j ) to (|A.1|) , the Polyakov phase probability distribution V(9) can be iden- 
tified with 



V(9) 



M 



1 d6 
1 



2tt 



1 



d9 



(A.13) 



in the approximation ( A.2j) . 

Assuming Z-i symmetry in the confined phase, we have from eq. ( 1.1U| ) 



V(9) 



2tt 



1 + 2 Y, P q cos(q9) 
g =even 



(A.14) 



which means that V(6) is even w.r.t. reflections around 6 = 0, ±|; i.e. 

V(tt-9) = V(9) 
V(--k-6) = P(9) 
V(-9) = V(9) 

Comparing eq. ( |A.15|) with (|A.13|) , we find that d~8 ph /d9 is also even under reflections 
arou] 
that 



(A.15) 



around 9 = 0, ±|. Since the derivative of an odd function is an even function, this means 



9 ph [9] = a + 



(A.16) 
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where 0(6*) is odd under reflections around = 0, and a is a constant. However, since 

9(6 M = 0) = 0^(0) = (A. 17) 

it follows that 6 ph [6 = 0] = f h (e M = 0) = 0. Then a = 0, and ^[^is odd around = 0. 
Further, from ( |A.8| ), (|A.9|) , and the assumed single- valuedness of 0(0 M ), it follows that 
M {6 = ±7r) = ±7r, and therefore that 

W h [±ir] = (A.18) 

Then, since 0^[O] = 9 pn [n] = 0, and d^/dO is even w.r.t. reflections around |, it follows 

that r [f] = 0, and that r [0] is odd w.r.t. reflections around |. By the same reasoning, 

W [0] is also odd w.r.t. reflections around — |. To summarize, ^ [0] has the reflection 
properties: 

^V^] = _r ft [0] 

P/l [7T-0] = -0 P "[0] 
TjP^ r Til TTP^ 



-7T 



0] = -0 pn [0] (A. 19) 



where the last two relationships are a consequence of global Z 2 symmetry in the confined 
phase. Therefore 

ef = - r /2 d9 M e p \e M ) 



IT JO 



4 r dQ V(Q)W h \B\ 
Jo 

4 dQ V{tt - 0)0^ [tt - 0] 

Jtt/2 

-4 r d~6 V(9)6 ph \S\ 

Jtt/2 

Jl r d e M r h (6 M ) 

IT Jtt/2 



= -efj (A.20) 

This explains why the correlation between 9 ph and 9 M found numerically in eq. ( |3.4| ) is a 
consequence of global Z 2 invariance. 

Our second task is to understand why Pi is negative in the 0' = + 9 M additive 
configuration, given that V{6) is peaked around = | as discussed in section 3. Introducing 
the probability distribution V'{9') for the Polyakov phases in the additive configurations 



< F(e') >= f w do' F{e')v\e') (A.21) 

J — 7T 
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and again neglecting the fluctuations of 9 ph at fixed 9 M around the mean 

1 ~/~ m i SP h laM\ 



<F>= r d9 M —F(29 M + 9 ph (9 M )) 



Under the change of variables 



0' = 26 M + T h (e M ) 



eq. (|A.22|) becomes 



(A.22) 
(A.23) 

, n , . (A-24) 

-2tt 27T d9 

The limits of integration have changed, but the original limits can be restored using the 
27r-periodicity of the integrand. To demonstrate the periodicity, we first have 

e\e + Ti) = e + n + e M (e + 7T) 



< f >-. 



2vr . 1 fjQM 

d9 ——z=rF{9 



9 + n + (9 + n -W h [9 + it\) 
9 + Ti + (9 + Ti + 9 ph [-9]) 



9 + 71 +{9 + 71 

9' '(B) + 2tt 



(A.25) 



where the reflection properties ( |A.19|) have been used. Single-valuedness of 9 (9) then 
implies the converse property 



9(9 + 2tt) = 9(9 ) + 7T 



Next, 



d9 



M 



dff 



d9 M d9 
1 d~9 



dr h " 1 



d9 



(A.26) 



(A.27) 



Then, applying ( |A.26| ) plus the fact that d9 ph /d9 is even w.r.t. reflections 9 — > —9 and 
9 -> ti - 0, 



'd9 hr 



-2n 



d9 



■ph ' 



d9 phK ^ 



d9 



d0 



d9 



dW h ^ 



)+7T 



d9 



' )+7T 



-1 



d9 )__, 



d9 

■ -9(f) 



(A.28) 
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Since F{9) is periodic by definition, this establishes the 27r-periodicity of the integrand in 
( |A.24j ), which can then be written 



< F > 



f 71 " -/ 1 dd M 




F{9 



(A.29) 



Comparing ( |OgD with (|A~21| ) 




de 



(A.30) 



Single-valuedness of 9 M (9) implies that d9 ph /d9 < 1, and with this restriction V{9) is 



a maximum where dW h /d9 is a minimum. However, we have previously deduced from the 
the fact that P\ = and P2 < that the probability distribution 

V{9) = 1 - 



dW h 



de 



(A.31) 



is Z 2 invariant and peaked at 9 = ±7r/2. Again, this distribution is maximized when 
d9 ph /d9 is minimized, which implies that dW h jd9 is a minimum at 9 = ±7r/2. Finally, 



9 9 = ±- 



71 



71 

2 

7T 



+ M ( 9 



2 V2 l 2 



7T 



(A.32) 



As a consequence, d9 ph /d9 is minimized at 9' = n (and also, by the same arguments, at 
9 — — 7r), which means that P(0 ) is peaked at 9 = ±tt, as illustrated in Fig. |15[ This 
explains why we expect Pi < in the additive configuration. 



References 

[1] G. 't Hooft, Nucl. Phys. B138 (1978) 1, 

G. Mack, in Recent Developments in Gauge Theories, edited by G. 't Hooft et al. 
(Plenum, New York, 1980) 

J. M. Cornwall, Nucl. Phys. B157 (1979) 392 

H. B. Nielsen and P. Olesen, Nucl. Phys. B160 (1979) 380; 
J. Ambj0rn and P. Olesen, Nucl. Phys. B170 (1980) 60; 265; 

J. Ambj0rn, B. Felsager, and P. Olesen, Nucl. Phys. B175 (1980) 349; 
R. P. Feynman, Nucl. Phys. B188 (1981) 479. 



35 



[2] 
[3] 

[4] 
[5] 

[6] 
[7] 
[8] 
[9] 
[10] 



G. 't Hooft, Nucl. Phys. B190 [FS3] (1981) 455. 

L. Del Debbio, M. Faber, J. Giedt, J. Greensite, and S. Olejmk, Phys. Rev. D58 (1998) 
094501, |hep-lat/9801027 . 



[11] 
[12] 

[13] 
[14] 
[15] 
[16] 
[17] 
[18] 



[19] 
[20] 
[21] 
[22] 



Ph. de Forcrand and M. D'Elia, Phys. Rev. Lett. 82 (1999) 4582, |hep-lat/9901020| . 

L. Del Debbio, M. Faber, J. Greensite, and S. Olejmk, in New Developments in Quan- 
tum Field Theory, ed. Poul Henrik Damgaard and Jerzy Jurkiewicz (Plenum Press, 



New York-London, 1998) 47, |hep-lat/9708023 



J. Ambj0rn and J. Greensite, JHEP (1998) 9805:004, |hep-lat/9804022 . 
A. Polyakov, Nucl. Phys. B120 (1977) 429. 
M. Zach, M. Faber, and P. Scala, Nucl. Phys. B529 (1998) 505, |hep-lat/97090T7 
A. Hart and M. Teper, |hep-lat/990203l| . 



M. Chernodub, S. Kato, N. Nakamura, M. Polikarpov, and T. Suzuki, [hep- 
lat/9902013; . 



L. Del Debbio, M. Faber, J. Giedt, J. Greensite, and S. Olejmk, Nucl. Phys. Proc. 
Suppl. 53 (1997) 141, |hep-lat/ 9607053 . 



H. Shiba and T. Suzuki, Phys. Lett. B333 (1994) 461, |hep-lat/94040T5| ; 

J. Stack, S. Nieman, and R. Wensley, Phys. Rev. D50 (1994) 3399, |hep-lat/9404UT4. 



T. De Grand and D. Toussaint, Phys. Rev. D22 (1980) 2473. 
P. Stephenson, Nucl. Phys. B550 (1999) 427, [hep-lat/9902002 



O. Philipsen and H. Wittig, Phys. Lett. B451 (1999) 146, |hep-lat/9902003 
C. Michael, Nucl. Phys. B Proc. Suppl. 26 (1992) 417. 



M. Faber, J. Greensite, and S. Olejmk, Phys. Rev. D57 (1998) 2603, |hep-lat/9710039 



M. Faber, J. Greensite, and S. Olejmk, Acta. Phys. Slov. 49 (1999) 177, [hep 1 
lat/9807008| . 



G. Poulis, Phys. Rev. D54 (1996) 6974, |hep-lat/96010T3 



G. Bali et al., Phys. Rev. D54 (1996) 2863, |hep-lat/9603012 . 
T. Suzuki et al., Phys. Lett. B347 (1995) 375, |hep-lat/9408003 



M. Chernodub, M. Polikarpov, and A. Veselov, JETP Lett. 69 (1999) 174, [hep- 
lat/9812012 . 



36 



[23] A. Hart and M. Teper, Phys. Rev. D58 (1998) 014504, |hep-lat/9712003] 

[24] L. Del Debbio, M. Faber, J. Greensite, and S. Olejmk, Phys. Rev. D55 (1997) 2298, 
|hep-lat/9610005 . 



[25] G. Bali, |hep-ph/980935ll . 

[26] J. Ambj0rn, D. Espriu, N. Sasakura, Mod.Phys.Lett. A12 (1997) 2665; For tsch. Phys. 
47 (1999) 287. 

[27] N. Evans, S.D.H. Hsu and M. Schwetz, Phys. Lett. B355 (1995) 475. N. Evans, 
S.D.H. Hsu, M. Schwetz and S.B. Selipsky, Nucl. Phys. B456 (1995) 205. N. Evans, 
S.D.H. Hsu and M. Schwetz, Nucl. Phys. B484 (1997) 124; Controlled Soft Breaking 
of N=l SQCD, |hep-th/ 9703 197 . 

[28] L. Alvarez- Gaume, J. Distler, C. Kounnas and M. Marino, Int. J. Mod. Phys. All 
(1996) 4745. L. Alvarez- Gaume and M. Marino, Int. J. Mod. Phys. A12 (1997) 975. 
L. Alvarez-Gaume, M. Marino and F. Zamora, Softly Broken N=2 QCD with Massive 
Quark HypermulUplets I, II, |hep-th/9703072] , 9707017. 

[29] O. Aharony, J. Sonnenschein, M.E. Peskin and S. Yankielowicz, Phys. Rev. D52 
(1995) 6157. E. D'Hoker, Y. Mimura and N. Sakai, Phys. Rev. D54 (1996) 7724. K. 
Konishi, Phys. Lett. B392 (1997) 101. K. Konishi and H. Terao, CP, Charge Frac- 
tionalizations and Low Energy Effective Actions in the SU(2) Seiberg-Witten Theories 
with Quarks, |hep-th/9707005| . 

[30] M. Engelhardt, K. Langfeld, H. Reinhardt, and O. Tennert, |hep-lat/ 9904004] ; and 
Phys. Lett. B452 (1999) 301, frep-lat/ 9805003 ; 

R. Bertie, M . Faber, J. Greensite, and S. Olejmk, JHEP 9903 (1999) 019, f^p- 
lat/9903023 . 



37 



